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^ '. Abstract 

' The integrals of motion of the tricritical Ising model are obtained by Thermodynamic Bethe 

Ansatz (TBA) equations derived from the A4 integrable lattice model. They are compared 



(N 



' with those given by the conformal field theory leading to a unique one-to-one lattice-conformal 



O 
(N 



O 



correspondence. They can also be followed along the renormalization group flows generated by 
the action of the boundary field fi^s on conformal boundary conditions in close analogy to the 
usual TBA description of energies. 



> . 1 Introduction 
O . 

The study of two-dimensional integrable quantum field theories (QFTs) and integrable lattice 
ly^ , models has lead to the discovery of a large number of deep analogies between them. It is worth 
I to recall here that, in both cases, integrability is equivalent to an algebraic relation (Yang-Baxter 
equation) for the fundamental "evolution" operator in the two languages, the S-matrix in the former 
case and the transfer matrix in the latter case. Furthermore, a possible link between them is given by 
the so-called conformally invariant field theories (CFTs) because they play a fundamental unifying 
Q^' role as being themselves integrable and describing the critical properties of lattice models. Hence, 
^ . adopting this point of view, one can consider integrable QFTs obtained by adding a perturbation 
which preserves the integrable structure of the conformal field theory. 

The relevant case, here, is the 99^ 3 boundary perturbation of the c = 7/10 unitary minimal 
^ ■ model (it describes the universality class of the tricritical Ising model (TIM)). Indeed, if this theory 
is defined on a cylinder of finite length, the operator (^1^3 acting on the boundaries generates a set 
of integrable flows [1] between different sectors of the model. 

On the other hand, many results indicate that the A4 Andrews-Baxter-Forrester (ABF) inte- 
grable lattice model [6] belongs to the same universality class of TIM. In particular, the boundary 
flows of TIM were extensively studied in [2, 3] using the continuum scaling limit of this model (see 
[4, 5] and references in [2] for the derivation of the TBA equations). Moreover, this correspondence 
has been made stronger by showing that the scaling limit of both the energy and the associated 
degeneracies of the lattice model exactly matches with the corresponding quantities given by the 
Virasoro characters of the c = 7/10 minimal unitary CFT [7]. 

The aim of this paper is to establish a lattice-conformal dictionary (the wording is borrowed 
from [8]) between the lattice model and the (possibly perturbed) underlying CFT. In other words we 
will find a correspondence between the states on the lattice, given by the eigenvalues/eigenvectors 
of the transfer matrix, and the states which belong to the Virasoro irreducible modules. In order to 
achieve this result we will look at a suitable family of involutive local integrals of motion (IM) . From 
the point of view of CFT such integrals of motion can be constructed as appropriate polynomials 
of the Virasoro generators whose actual expression can be obtained in a constructive way (since 
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a general expression is not known to exist), sec [11]. Within the lattice model, such integrals of 
motion can be accessed by usual TBA techniques showing a remarkable correspondence between 
their eigenstates and those of the transfer matrix (a similar analysis was done in [12] by means of 
the nonlinear integral equations method in the context of KdV theory.). Therefore, as we will show 
in the sequel, the desired statc-by-state correspondence can be found by comparing the eigenvalues 
of the integrals of motion in the CFT with those of the lattice model. In this sense, our approach 
is algebraic-analytic while in [13] it is combinatorial. 

Finally, one can wonder about the fate of such integrals of motion when a relevant perturbation 
is switched on. In the present case, where the (pi^^ boundary perturbation is concerned, the answer 
is that they can be chosen to be preserved, i.e. to remain involutivc integrals of motion under 
such a perturbation. Indeed, among the possible families of involutive integrals of motion allowed 
in the CFT, we will consider that one whose spins are predicted to be preserved, for the present 
perturbation, by the counting argument in [14]. This implies that the perturbed theory is still 
integrable and it is remarkable that the eigenvalues of the integrals of motion can be followed along 
the Renormalization Group flow in close analogy to the usual TBA description of energies. 

As a final consideration, we also find a connection between our results and a recent series 
of papers by Bazhanov, Lukyanov and Zamolodchikov (BLZ) [16, 17]. We will see that, up to 
normalizations, the eigenvalues of the transfer matrix (and hence those of the integrals of motion) 
of the lattice model coincides with the eigenvalues of the so-called "T-operators" introduced in the 
above mentioned papers. 

The outline of the paper is as follows. In sect. 2 we describe the A4 ABF integrable lattice 
model; in sect. 3 we carefully discuss both the continuum scaling limit of the lattice transfer 
matrix and the emergence of the integrals of motion in such a context; sect. 4 is devoted to a brief 
presentation of the integrals of motion in CFT. We give our results in sect. 5. The connection 
between our paper and the BLZ papers is elucidated in sect. 6, and finally we draw our conclusions 
in sect. 7. 



2 A4 lattice model 

This model was exhaustively discussed in [7, 2] and we summarize here the properties that will 

be useful later. For convenience, we will adopt the notations introduced in [2]. At the critical 
temperature, on a cylinder of finite length, it is defined by the following double row transfer matrix 



D^KOa,<T' = Di^l|2,l(^>6att)a,<7' = ? 2 
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(2.1) 



which has been specialized to the boundary flow that we are going to examine: xi,2 '-^ Xi,i (note 
that it is a constant r flow). 

The left boundary is kept fixed and the boundary interaction, shown with dashed lines, acts 
on the right boundary only and it is tuned by the parameter ^. The heights aj,aj,Tj G {1, 2, 3, 4} 
on adjacent vertexes are constrained to follow the A4 adjacency rule which also forces the number 
N of faces on a row to be odd. In the following we will always restrict to this case, in particular 
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when the hmit limjv-»oo will be considered. The bulk weights are fixed to their critical values 
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where the physical range of the spectral parameter is < n < A, and A = ^ is the crossing 
parameter. The crossing factors are defined as Sa = sinaA/ sin A. 

Among the general expressions for integrable boundaries given in [9], we only need to consider the 
following case: 
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The second line of the equation shows that the same boundary interaction can be obtained by 
adding a single column to the right of the lattice; the faces of this column depend upon ^ but the 
new right boundary is fixed to the sequence 1, 2, 1, . . .. 

Let us now consider the right boundary of the lattice defined in (2.1). We can see that, for 
^ = —GO it is formed by the quasi-free sequence 2, 2 ± 1, 2, 2 it 1, . . ., while for ^ = it is formed 

by the fixed sequence 1,2,1,2, In the continuum scaling limit, the former case corresponds to 

the ultraviolet (UV) critical point Xi,2{q) and the latter to the infrared (IR) critical point Xi,i(q') 
of the boundary renormalization group flow. 

By construction, the given bulk and boundary weights satisfy the Yang-Baxter and Boundary Yang- 
Baxter equations leading to an integrable model. Consequently, the double row transfer matrix 
forms a one parameter family of commuting operators: 



D^(t.,O,D^(u',O]=0 



(2.4) 



for arbitrary complex values of the spectral parameters u,u'. An important symmetry becomes 
apparent when one defines the following normalized transfer matrix 



sin(ii + 2A) sin A 
_sin(« + 3A) sin(M + A)_ 



1 2N 



where 
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5'2(u,6att) 
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sin(2ti + A) sin(2ii — 3A) ' 

- sin A s'm{u - ^latt + A) sin(u + ^latt + 3A) cosh 2Im(^iatt) 

sinfn 



6att) sm{u - ^latt + 2A) sin(u + ^latt - A) sm{u + ^latt + 2A) 
Indeed, one can show that the following functional equation is satisfied: 

t^(n, t^(tx + A, = 1 + t'^iu + 3A, 0- 



(2.5) 

(2.6) 
(2.7) 

(2.8) 
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Fig. 1. Schematic example of zeros of the eigenvalue of the transfer matrix labeled (2 1 1 0|1 0). Here: 

mi = 6, 7712 = 2, m = 3, 712 = 1- 



The normalized transfer matrix t also forms a family of commuting operators 

[t^KO,t^(^',0]=0 (2.9) 

so that the eigcnstates are independent of u, and the functional equation holds true also for the 
eigenvalues. We would like to stress that it is a crucial point since the TBA approach is precisely 
based on such an evidence. 

In the continuum scaling limit, the large N corrections to the eigenvalues^ of the double row 
transfer matrix are related to the excitation energies of the associated perturbed conformal field 
theory by 

-ilogZ^^KC-logiV) = Nh{u) + fUu,0 + ^^^^^m)+oi^) (2.10) 

where is the bulk free energy, /^^ is the surface (i.e. boundary dependent) free energy and E{^) 
is a scaling function given by (3.20). Since both the bulk and surface free energy contributions are 
the same for all the eigenvalues, they can be removed from the previous expression by subtracting 
the largest eigenvalue of D'^ (it plays the role of a hamiltonian ground state). At the boundary 
critical points, the scaling energy E{^) reduces to 

m) = \ , , (2.11) 

where n, n' G N are certain excitation levels. 

The derivation of both the TBA equations and the scaling energy (2.11) is obtained by looking 
at the analytic structure of the eigenvalues of the transfer matrix: D^{u,^) is an entire function 

of u characterized by simple zeros only, moreover it is periodic on the real axis u = u + n. In the 
large N limit, the zeros arc organized in two sort of strings and they are distributed in two different 
strips, defined by — ^ Re{u) ^ |A and 2A ^ Re(u) ^ 4A respectively. An example of such a 
situation is schematically shown in Fig. 1. 

The zeros are bound either to appear in the center of each strip (located at Re(u) = ^ or 3A) and 
we will refer to them as 1-strings, or they can appear in pairs {u, u') with the same imaginary part 
Im(u) = Im(ti') on the two edges of a strip, and we will call them 2-strings. In order to label the 
number of such strings in the upper-half plane, we will use mi and m-2 for 1-strings located in 

^We use D, t to indicate eigenvalues of D, t. 
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the first and second strip respectively; analogously ni and n2 will be used for 2-strings^. These 
numbers form the so-called (m, n)-system and are constrained to be 

m = - mi ^ 

712 = ^ - m2 ^ 

The imaginary part of the position of the 1-string is indicated by v^' 

first strip: {^X + iv'j^\^) = 0, second strip: {3\ + iv''^\ ^) = 0; (2.13) 

(7) (7) 

where the /c's are ordered from the bottom to the top in the upper half plane, < v)^ < ■Ug < 
. . . < Vm\ . Analogously, w^^^^ is used for the 2-strings. 

Finally, the states (cigcnvalucs-cigcnstatcs) both on the lattice and in the continuum scaling 

(7) 

limit are uniquely characterized by the non-negative quantum numbers Ij: . Their topological 
meaning is that when we refer to the k-th 1-string in strip j, I^' is the number of 2-strings with 
larger imaginary part. The notation which indicates the content of zeros works as in the following 
example as shown in Fig. 1: (211000|10) is a state with 6 zeros in the first strip and 2 in the 
second; their quantum numbers are respectively = 2, ll^^^ = 1,.... In the sector (1,2) a 
frozen zero can occur [7], so a parity cr = ±1 is added to the previous notation (namely ( • | • )±, see 
Tables 4 and 5). 



mi,m2 even. (2-12) 



3 Continuous transfer matrix from the lattice 

Since our task is to compare lattice computations with CFT results, we have to define a suitable 
continuum limit of lattice quantities. In the following we will first define the continuum scaling limit 
for the lattice transfer matrix, and then the corresponding integrals of motion will be introduced. 



3.1 Continuum scaling limit of the transfer matrix 

In [2] the finite size (order 1/N) corrections to the lattice energies (scaling energies) were computed 
with standard TBA techniques that involve a continuum scaling limit on the lattice model. This 
becomes fully apparent when the bulk of the lattice is off the critical temperature, t ^ {t = 
(T — Tc)/Tc), and the continuum limit (N 00, a ^ 0, a is the lattice spacing) is sensible only 
if the temperature is also scaled, t as N\t\'^ = = const, {y = 5/4) [10]. Physically, this 
corresponds to a huge magnification of the critical region. The bulk critical point itself is reached 
when the dimensionless mass ji goes to zero: ^ In the present case, the bulk is at the critical 
temperature (/x = 0) and we use the boundary perturbation to move away from the boundary 
critical points. 

The aim of this section is to show how to define a meaningful continuum scaling limit of the 
lattice transfer matrix. Let us recall that in [2], from the derivation of the TBA equations, it was 
shown that the following scaling property for the normalized transfer matrix holds 

t(u, = lim t^(u + l\ogN,i- log A^) (3.1) 

and hence a proper definition of the continuous transfer matrix (i.e., a transfer matrix for the 
continuum theory) is 

2N 

(3.2) 



D(n, = lim D'^ (n + - log A^, ^ - log N) 

AT— »oo 5 



2 sin Ae*("-t) 



^We do not need to keep track of the lower half plane because it is the mirror image of the upper half plane, 
thanks to the complex conjugation symmetry of the transfer matrix. 
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where we can explicitly observe the periodicity property D(u, ^) = D(u + tt, ^). Since is real 
in the center of each strip (as shown in [2]) and the square of the factor in brackets is real for 
Re(?x) = A/2, 3A, we deduce that also D is real in the center of each strip. 

The main motivation to introduce the previous definition is to have an object for the continuum 
theory with the same analytic structure of the lattice transfer matrix. Hence, let us analyze such a 
structure starting from the string content of the eigenvalues of the double-row transfer matrix. A 
given eigenvalue D^{u, ^) is characterized by its content of 1-strings. It is important to notice that, 
increasing N, the number of 1-strings, of 2-strings in the second strip and the relative position of 
1- and 2-strings docs not change. Instead, it is observed the appearance of new 2-strings close to 
the real axis in the first strip, moreover eq. (2.12) tells us that ni grows as N/2 (notice that N 
grows in steps of 2 so that ni grows in steps of 1). Such new 2-strings push away the remaining 
zeros from the real axis and since the extension of the region occupied by them grows as logA'^, 
then the following limit is finite {v^^ is a function of N): 



(j) def ^r- (j) 



logiV). 



(3.3) 



In terms of the scaled coordinate , the real axis located at Im(tt) = is shifted to —ioo. From 
(3.2) and (3.3) one is immediately led to the following observation: 



hm i?^(lA + ^(54^) 

N-^oo 2 5 



log AT) + - log AT, ^ 



logN) = Q = D{\\+'-y\^\i) (3.4) 



(where a similar equation holds for the zeros in the second strip). The remarkable consequence of 
such an equation is that the content of zeros of an eigenvalue survives the continuum limit, the 
only difference being the infinite number of 2-strings that appear in the first strip, for u — ^ —ioo, 
because of the growth of rai ~ A^/2. 

In order to show that the limit in (3.2) is indeed meaningful, in the sense that it is finite, we have 
to take into account the following ingredients. Firstly, the derivation of the TBA equations shows 
the existence of the limit in (3.1); secondly, examining the normalization coefficients S'2('u,^), S{u) 
in eq. (2.5) we have 



S{u) 



lim S2(n + J log AT, -A + (e - log N)) 

N-^oo 5 5 



- sinA sm{u + f + 2A) 
sm{u + - 2A) sin(u + + A) ' 

lim ^(M-h -logA^) = 1 



(3.5) 



(3.6) 



which are finite quantities. The factor in square brackets in (2.5) gives: 

27V 



sin(-;i + ^logN + 2A) sinA 



2sinAe^("-2) 
iVV5 



2N 



sm{u I log A -h 3A) sm{u | log A -|- A) 
that exactly cancels the corresponding factor in (3.2). Therefore we get the equation 

t(^x,0=D(n,0^2(n,0 



(3.7) 



(3.8) 



which shows that D(u, ^) is finite because i{u,$,) and S2 are finite. It also shows that, in the 
continuum theory, t and D are almost equivalent, the difference being in zeros and poles explicitly 
dependent upon ^ that are introduced in t by the factor ^2. On the other hand, the lattice 
possesses some non physical zeros and poles of order 2N contained in the square bracket term of 
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(2.5). Moreover, for a generic Ap ABF model, the zeros have to be understood as zeros of 
while their occurrence in is much more involved. 

Finally, it can be shown that the functional equation (2.8) keeps the same form even after the 
continuum scaling limit 

i{u, i{u + A, = 1 + 1(« + 3A, 0- (3.9) 

This equation can be solved for the eigenvalues of t and the solution is given by the TBA equations 
obtained in [2]. 

It is important to notice that the normalization of t^{u, ^) is uniquely fixed by the requirement 
that the functional equation holds in the form (2.8), namely no multiplicative or additive constants 
are allowed (in particular, this fixes the relation (2.5) between D^('u,^) and t^(u, ^)). Since such a 
functional equation remains unchanged after the continuum scaling limit, the lattice normalization 
uniquely fixes the continuum one. This fact will have important consequences for a correct definition 
of the integrals of motion. 

3.2 TBA equations and integrals of motion 

Let us discuss the TBA equations for the eigenvalues of t. First of all, we define the functions 
ti{x,^) and t2{x,$,) which are the eigenvalues of t referred to the center of each strip: 

= ii^+'-x,0 = s,e-^^^^), (3.10) 

t2{x,0 = i(3A + ^x,0 = 526-^2(2;). (3.11) 
5 

where the pseudoenergies ej{x) are introduced, and the quantities Sj = ±1 play the role of inte- 
gration constants and will be fixed later. Since the functions ij{x) are real for real x and satisfy 
tj{x) > —1, we can define the real functions Lj{x) as 

Lj{x)=log{l + ij{x,0)- (3.12) 
In terms of the functions ti and t2, the functional equation (3.9) becomes 

ti{x + i^) ti{x - i^) = l + i2{x) 

t2(x + i|)i2(x-i|) = l + ii{x). (3.13) 
Consequently the TBA equations can be written as 

h(.,() ^ gta^h^ exp (/;;|_^) , (3.14) 

= e--%.&(..9nta^^^exp(/;;|-^) (3.15) 

and the positions of the zeros are fixed by the quantization conditions (note the inversion of the 
indices: ipi is for strip 2 and '02 for strip 1): 

MVk^) = 4^^^ = (2(4^^ + mi - A;) + 1 - ma) tt, (3.16) 

MVk^) = 4^^7r = (2(/f ) + m2 - fc) + 1 - mi) tt. (3.17) 
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The same equations hold for the 2-string locations so, each time ipj (x) = mr is satisfied for an 
integer n with the appropriate parity^, a; is a 1-string or a 2-string in the strip 3—j. The explicit 
expressions for i/jj are given by: 

^i(x) = Re(-iei(x-i|)) =Re(zlog(siti(x-z|,0)) (3.18) 

= i\oggAx-i-,i)+iy logtanhf^^- \- t-r) - I tt- • , / 

2'^^ ^ ^ ^2 4^ W-oo 27r sinh(x - y) ' 

ip2{x) = Re ^-ie2(x - i^)^ = Re ^ilog(s2t2(a; - (3-19) 



4e-^ + zlog^2(a^-«|,0 +^Xl^°Stanh(^^— ^ ^^^/ 

k=l ^^J-i 



27r sinh(x — y) 



where the integral around the singularity x = y is understood as a principal value. Moreover, 
for numerical convenience, we take the fundamental branch for each logarithm so that in general 
logo + log 6 cannot be identified with log(a6). The energy predicted for each state is given by: 

m = - E - / 4 (3.20) 

and reduces to (2.11) at both the critical points (compare also with (2.10)). Since we are interested 
in the flow xi,2 ^ Xi,i) we have 

si = S2 = 1; 

X + £ (3.21) 
5i(a;,0 = l; 52(a;,0 =tanh^— . 

It is useful to recall that the functions ij{x, £) are eigenvalues of a commuting family of operators: 

[i{u,O,i{u,O]=0. (3.22) 

Hence we are in the position to show how the local integrals of motion emerge from the continuum 
scaling limit of the transfer matrix t(u,^). First of all we expand the quantity logti for x — oo 
where ti is given by (3.14). In order to do so we need the following formulae 

yi}^ — X 1 — e^^fe 
log tanh -5^- = log — 

1-t / 



log^ = -2(t + - + - + ...), |t|<l 

oo 

^ = 2y(-l)"-ie(2'*-i)^, x<0. 



cosh X 

n=l 

Hence, taking advantage of the previous results, we arrive to the following expressions for the 
continuum scaling limit of the logarithm of the transfer matrix 

oo 

logti(a;,0 = -j;C„/2n-i(0e^'"~'^", (3.23) 

n=l 

Cnhn-m = ;^f:e-(^"-^)^i%(-l)" r"^e-(2"-i)^L2(y), (3.24) 



^There is the following constraint: v}^' + — = mod 2 
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where the quantities^ hn-i arc the cigenvahies of the infinite family of integrals of motion in 
involution {/2n-i}n- It is important to stress that the TBA ealeulation provides only the value of 
the product Cnhn-iiO ^ is clear from eq. (3.23) and (3.24). The coefficients C„ will be fixed 
later by imposing appropriate normalization conditions (see sect. 5). 

As pointed out in sect. 3.1, the normalization of ti is fixed once forever on the lattice and 
survives the continuum scaling limit hence its largest eigenvalue, being the lattice groundstate, in 
that limit is directly mapped to the true vacuum of the continuum field theory. This forbids the 
presence of additive constants to hn-i in (3.23). Such a thing is in agreement with the general 
properties one expects for the above integrals of motion in the continuum limit: since they are 
given by integrals over local densities, and the local densities are components of Lorentz vectors, 
there is no freedom to add any arbitrary scalar contribution to them. Hence we expect to find, at 
the fixed points of the massless flow, the same values for hn-i given by the conformal field theory 
(see next section) once the proper normalization for the corresponding C„ has been chosen. 

As an example, it is interesting to note that comparing (3.24) with (3.20) we obtain Ci = tt 
and hiO = ^iO- expected, the first integral of motion is the energy. 

Now, our purpose is to use the TBA equations to compute the eigenvalues of /2n-i and 
compare them with their analog from the conformal field theory. As we will see, the two sets of 
integrals of motion precisely correspond to one another. 



4 Integrals of motion in conformal field theory 

Conformal field theories are integrable, in the sense that they possesses an infinite number of 
mutually commuting integrals of motion. Their construction was given in [11] by canonically 
quantizing the family of classical integrals of motion of the Modified Korteweg-de Vries equation 
(that also coincide with the integrals of motion of the classical sine-Gordon equation). 

A general expression for such integrals of motion is not known but they can all be obtained 
in a constructive way. They are local expressions (polynomials) in the Virasoro generators and the 
first few of them are given by^ (see [16] for their explicit expressions) 

Ii = Lo-^ 



I 



3 



24 

2±L.„L„ + Ll-^-±^L„ + '-^^^ (4.1) 



12 2880 

n=l 

3 °° 

m.n,p n=l 

°^^11 + C9 c \^ ^ 4 + c, 



n=l 

(2 + c)(20 + 3c) _ c(3c+ 14)(7c + 68) 
^ 576 ° 290304 

(the symbol : : indicates normal order, obtained by arranging the operators L„ in an increasing 

sequence with respect to their indices). 

Their simultaneous diagonalization is, in principle, a straightforward but lengthy calculation 
of linear algebra. Actually, the rapidly growing complexity of the involved expressions makes the 

''We hope that the symbols l''f^\ with j = 1,2 introduced for the quantum numbers, arc not confused with the 
similar symbols /2n-i used for the integrals of motion. 

®It is convenient to indicate each integral of motion with /2n-i and the corresponding eigenvalues with hn-i- 
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computation possible for the first few levels only. This becomes clear if we look at the characters 
of the Virasoro algebra: 

where the trace is taken in the irreducible module of the highest weight state (r, s) . For example, 
for the TIM one has: 

XiAl) = 9"^(l + g^ + g^ + 2g^ + 2g5 + 4^6 + V + 7gS + ...) (4.2) 
Xi,2(g) = + g + + 2^3 + 3^^ + 4g5 + 6g6 + 8g^ + ...). (4.3) 

In the previous expression, the coefficient in the monomial anQ^ gives the degeneracy of the 
subspace at level n. At a given level of truncation, the expression for the vacuum sector Xi,i('?) 
holds for all the unitary minimal models except the c = 1/2 Ising model. Up to level 3 there is no 
degeneracy and the eigenstate is given by the unique linearly independent vector. At level 4 and 5 
the double degeneracy leads to an algebraic quadratic equation that can be easily solved. At level 
6 the degeneracy of order 4 leads to a quartic algebraic equation and, in general, one realizes that 
the complexity in managing this problem rapidly grows with the degeneracy of the level. More 
details on this topic can be found in Tables Al and A2 of the Appendix where we present a list 
of eigenstates/eigenvalucs for the vacuum sector of all unitary minimal models (Table Al) and for 
the sector (1, 2) of TIM (Table A2). From Tables Al and A2, we notice that alone is enough to 
completely remove the degeneracy, at least up to level 6. 

The list of integrals of motion discussed here survives the off critical perturbation by (pi^s [14]. 
More precisely, off critical involutive integrals of motion /2ri-i(A) exist for the perturbed CFT^ 
Mp^q + Xipi,3 and at criticality they reduce to the expression (4.1). 



5 TBA data versus CFT 

As anticipated in the introduction, we are interested in studying the model both at criticality 
and along the boundary flow xi,2 ^ Xi,i generated by the (/?i^3 integrable perturbation of TIM. 
Let us first describe the sector (1,1) of the theory, i.e. the IR critical point. Let us recall that 
TBA computations provide the product Cnhn-iiO^ ^l- (3.24). The constants C„ are precisely 
introduced to make contact with the integrals of motion of the CFT. As already pointed out at the 
end of sect. 3.2, the first of them is exactly know, Ci = n. Since we are interested in J3 and J5 we 
can numerically fix the corresponding constants C2 and C3 using the vacuum state^ 

C2 = 2.1838434, C3 = 3.7555032. (5.1) 

The general expression was computed in [16] and reads^: 

Cn = 2^-" 3^-2" 5I-" TT. (5.2) 

" n!(4n-2)! ^ ^ 

In particular, this gives: Ci = tt, (72 = vr ~ 2.1838432 and C3 = |i§g§ vr - 3.7555026. 
After that, we are ready to turn the crank of TBA equations in order to obtain the other eigenval- 
ues of the operators 1 3 and J5 . The comparison between them and the corresponding eigenvalues 
coming from CFT can be found in Table 1. We can immediately notice that the agreement is quite 
remarkable: the solutions of the TBA equations are given with seven significant digits of precision 
and perfectly match the CFT results. 



^With Mp^q we refer to the minimal CFT with central charge c = 1 — ^^''^^^ ■ 
Actually, one reaches higher numerical precision using the first excited state and so we did. 
®We would like to thanks the referee of Nucl. Phys. B for having pointed out to us this equation. 
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Table 1. Comparison between the eigenvalues of J3, J5 from CFT and from TBA in the vacuum sector 
of TIM. The left column indicates the level degeneracy (l-d.) predicted by the character of the Virasoro 
algebra. 



1 A 

i.d. 


lattice-conformal dictionary 


J3 [Lb i j 


I3 (iBA) 


is (CFT) 


I5 (iBA) 


1 


- 


- |o) 


0061979 


0061979 


-0 0028301 


-0 0028301 




(00) ^ 


L-2\0) 


4.2561979 


4.2561979 


8.1731074 


8.1731074 


1 ^3 
Lq 


(10) <- 




ly.ioiiyo 




1U4. 14040 


1 r\A 1 A^A^ 
1U4. 14040 


2/ 


(20) ^ 


3(^±^L_4 + 2L2 2)|0) 


52.052045 


52.052042 


556.20159 


556.20155 




(11) ^ 


3(4-f^L_4 + 2L2_2)|0) 


22.560351 


22.560348 


108.29650 


108.29648 




(30) ^ 


(I±^^L_5 + 20L_3L_2)|0) 


110.13688 


110.13687 


1943.8246 


1943.8244 




(21) ^ 


_^ (I=^^^_^^20L_3L_2)|0) 


55.125517 


55.125511 


551.34946 


551.34937 




(40) ^ 


(11.124748 L_6 + 9.6451291 L_4L_2 
+4.4320186 L?.3 +L3 2)10) 


200.49775 


200.49773 


5280.2710 


5280.2705 




(31) ^ 


(-4.9655743 L_6 + 2.3354391 L_4L_2 
+0.71473858 L?.3 +L_2)|0) 


112.78147 


112.78146 


1919.0972 


1919.0969 




(22) ^ 


(0.66457527 L_6 - 1.2909210 L_4-L_2 
-1.2605013Li3 + Li2)|0) 


69.265146 


69.265138 


636.39114 


636.39101 


(0000|00) ^ 


( 1.6612491 L_6 4.0646472 L_4L_2 
+1.4118691 L2 3 + Li2)|0) 


35.980431 


35.980428 


209.28433 


209.28430 
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Table 2. Comparison between the eigenvalues of I^, J5 from CFT and from TBA in the (1, 2) sector of TIM. 
The left column indicates the level degeneracy (l.d.) predicted by the character of the Virasoro algebra. 



l.d. 


lattice-conformal dictionary 




/3 (CFT) 


h (TBA) 


h (CFT) 


h (TBA) 


1 


(0)+ - 






-0.0063021 


-0.0063020 


0.0026543 


0.0026543 




(0)_ ^ 


- ^-ilio) 




1 3fi8fiQ7Q 


1 3686970 


1 3062480 


1 3062466 


1 2 


(1)- - 






10.243698 


10.243694 


37.203592 


37.203567 




(2)_ ^ 


(l6L_3 + (3 + V649)L_2L_i) 


1 10' 


33.732663 


33.732653 


271.05745 


271.05732 




(000)+ <- 


(l6L_3 + (3- V649)i^-2i^-i) 


ll^) 


10.804733 


10.804732 


34.613169 


34.613163 




(3)_ ^ 


(L_4 + 2.1916731 L_3L_i 
+1.2645819 L?.2)| 3^) 




78.939940 


78.939920 


1117.9764 


1117.9759 




(100)+ ^ 


(L_4 - 0.50977440 L_3L_i 
-0.42935709 L2 2)13^) 




33.642011 


33.642008 


255.64791 


255.64786 


(000|00)_ ^ 


(L_4 + 2.5725923 L_3L_i 
-1.7645661L2_2)|X) 




16.699143 


16.699141 


64.926784 


64.926773 




(4)_ ^ 


(L_5 + 1.6541786 L_4L_i 

+1.3195562 L_3L_2 + 0.52470182 Lia^^-i 


152.96762 


152.96759 


3367.1124 


3367.1113 




(200)+ ^ 


(L_5 - 0.29373082 L_4L_i 
-0.012595542 L_3L_2 - 0.45966766^2 2 L 


78.030007 


78.029998 


1071.0704 


1071.0702 




(110)+ ^ 


(L_5 + 0.19041567 L_4L_i 

-2.9757621 L_3L_2 + 1.8161567^2 2^^-1) 


41.633949 

h) 


41.633947 


283.29346 


283.29343 


(100|00)_ ^ 


(L_5 + 4.9698982 L_4L_i 

-1.1753145 L_3L_2 - 1.4399411 L2 2L_i) 


43.143217 
1^) 


43.143213 


365.24996 


365.24990 



12 



Table 3. Flow xi,i Xi,2 (reverse of the physical flow). We present the explicit mapping of states from IR 
to UV. Here n}^, n^^ are the excitation levels above the ground states, respectively A = and A = 1/10. 



^IR 


Mapping of states 




^IR 


Mapping of states 


^IR 





- (0)+ 





5 


(21) ^ (100)+ 


4 


2 


(0 0) ^ (0)_ 


1 


6 


(40) ^ (4)_ 


5 


3 


(10) ^ (1)_ 


2 


6 


(31) 1-^ (2 0)+ 


5 


4 


(2 0) ^ (2)_ 


3 


6 


(2 2) ^ (110)+ 


5 


4 


(11) ^ (0 0)+ 


3 


6 


(0000|00) ^ (000|00)_ 


4 


5 


(30) ^ (3)_ 


4 









Hence we can draw the following conclusions: firstly, the expansion postulated in (3.23) is the 
suitable one for the computation of the eigenvalues of the integrals of motion in the lattice model 
(in the continuum limit) and shed some light upon the connection between the integrable structure 
of the CFT and the lattice model; secondly, such a matching between the eigenvalues allows to 
conjecture a one-to-one correspondence between the states of the CFT and the corresponding states 
defined within the lattice model. 

Few remarks are in order with respect to the latter issue. Even if we are able to build the cor- 
respondence of states level by level, we are not able to write down explicit, closed expressions for 

the eigenvalues /2n-i at criticality as it was done for the energy (sec (2.11) and more complete 
expressions in [2]). Notwithstanding this, numerical evidence strongly supports the conjecture that 
the state-by-state correspondence is exact giving a true lattice- conformal dictionary. 

Let us study the (1,2) sector of TIM, as in Table 2. As expected we get a similar result as 
before, confirming the above conjectured correspondence (it is useful to stress that the constants 
Cn are fixed once forever, because they are computed on the true vacuum that is independent of 

0- 

Finally, we can explore by means of TBA equations what happens to the integrals of motion 

when the boundary 991.3 perturbation is switched on (see [15] for the case of a massive perturbation 
in the context of [12]). Since off-critical expressions for the integrals of motion are not available we 
cannot make any comparison, but we can just follow the evolution along the whole boundary flow 
Xi,2 1-^ Xi,i- 

In Fig. 2 the first few eigenvalues of 1 3 are plotted (against ^): we can notice the absence of 
degeneracy and the clear correspondence between the states belonging to the UV and IR fixed 
points. Such a correspondence has been explicitly shown in Table 3. 

The method to solve numerically the TBA equations was widely discussed in [2] . The compu- 
tation of the integrals of motion does not require any new technique: it simply uses the previously 
obtained solution in the expressions (3.24). 



13 



100 - 



80 



60 - 



40 



20 - 



(1,2) 




(1,1) 



Fig. 2. Behaviour of the eigenvalues of the integral of motion /3(^) along the renormalization group flow. 
We point out that, in the UV region, there is no degeneracy of states; what seems a degeneracy is a graphical 
effect because of the vicinity of the eigenvalues (compare with Section 5). 
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6 Continuous transfer matrix from CFT 



In the present section we will discuss the relation between the so called T-operators introduced 
by Bazhanov, Lukyanov and Zamolodchikov (BLZ) in [16], and the continuum limit of the lattice 
transfer matrix t^^-* (3.1). 

In the series of papers [16], the authors took advantage of the integrable structure of Conformal 
Field Theories building up the quantum version of the monodromy matrices of KdV theory. 

A crucial role in their construction is played by the " T-operators" , which are the trace of such 
quantum monodromy matrices. As shown by BLZ, as a consequence of the Quantum Yang-Baxter 
equations the T-operators form an infinite set of operators in involution 

[Tj{u),Tjifi)]=0 (6.1) 

where u, /x are spectral parameters and the indexes j, j' label the corresponding representations 
of Uq{sl{2)). Such a quantum group structure has its root in the Feigin-Fuchs construction related 
to the quantum version of the monodromy matrices. Here it is important to recall that since the 
stress energy tensor is defined as 

-P^T{u) =: ^'{uf : +(1 - P^)ct^" {u) + ^ (6.2) 

where (j) is the compactified bosonic field (see [16] for details), then the relation between the central 
charge c of the CFT and the parameter (5 is 

c = 13 - 6(^^ + ^~^) (6.3) 

and the deformation parameter q of the quantum algebra is related to (5 as 

q = e^^^'. (6.4) 

BLZ analyzed in detail the properties of the operators T^jiy) establishing, in particular, that they 
satisfy the following (infinite) set of functional relations 

T,{q^u)T,{q--2y) = I + T .^,{v)T (6.5) 

showing that at particular values of the parameter (3 a truncation of such a system of functional 
relations occurs. In [16] this phenomenon was shown explicitly for the non-unitary series of minimal 
models M.2,2n+Z, giving rise to a set of equations which is equivalent to the F-system of the TBA 
equations for the massless S'-matrix associated to the previous CFT's. As pointed out in [16], 
the functional relations (6.5) also hold in the case of the (^1^3 perturbation of A^2,2n+3- Taking 
advantage of this fact in [17] the solutions of the TBA equations for the excited states of the Yang- 
Lee model (A^2,5) with (^1^3 perturbation were found. 

As anticipated in [16], such an analysis can be repeated mutatis mutandis for the series of minimal 
unitary models A1p^p+i where p = 3, 4, . . . is an integer which labels their central charges 

(6.6) 



p(jp+iy 



It can be shown that in this case the truncation restricts the range of variation of the index j of the 
Tj(i/) operators to the subset j G {0, ^, 1, . . . , leading also to the following symmetry relation 

T^ .(i/) = T,(i/). (6.7) 

2 
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If we refer to Tj {u) as the eigenvalues^ of Tj (u) , then the following relations can be written 

Tj{q'^u)Tj{q--2u) = I + T^^{iy)T^_i{iy) 
Te^-M = Tj{v); Tj^{iy) = To{u) = 1. (6.8) 

2 J 2 

Inspired by the results of [16, 18, 19], we perform the substitution (a = 2j) 

Ta(9) = Ta,i(iy)Ta_i(u); i/ = e^'^ (6.9) 

2"'"2 2 2 

which allows to cast eq. (6.8) in the form 

T„(e + ^)T,(e-^) = (l + T,+i(^))(l + T,_i(^)) 

Tp_i_„(^) = Taiey, Tp_,{e) = To{e) = o. (6.10) 

where a G {0, 1, 2, . . . , p — 1}. Now we are in the position to make contact with the system of TBA 
functional equations for the 1^91.3 perturbation of the theories Mp^p-^-i, obtained from a standard 
scattering approach [20]. Such equations can be written as 

p-i 

Yai0 + '-^)Yaie - f ) = 11(1 + n(e))^- = (1 + y„+l(^))(l + Ya-^ie)) (6.11) 

&=0 

where iab = 5o+i,b + ^a-i,b is the incidence matrix of the Dynkin diagram of Ap-2 and Yq(9) = 
Yp-i{9) = 0. The equivalence between (6.10) and (6.11) can be immediately recognized. 

The same equations were first obtained by Kliimper and Pearce [19] for the integrable RSOS 
fusion hierarchies. This reveals that there are identical structures coming from the integrable lattice 
and from the minimal models of conformal field theory. In both the cases a transfer matrix can 
be defined such that it obeys the same functional equations. The connection between the BLZ 
formalism and the transfer matrix formalism of [19] is given by the conjectured correspondence: 

Ta(e) ^ ta (-^^) (6-12) 

which means that they are equal modulo some lattice-continuum normalization constants. Such a 
conjecture has been numerically tested (at critieality) in sect. 5 for the tricritieal Ising model. In 
particular, we checked that the two terms in (6.12) contain the same integrals of motion eigenvalues, 
under the expansion defined in (3.23). 

The previous reasoning suggests that a suitable prescription to find the solutions to eqs. (6.10) 
can be that of determine the analytic structure of the by means of the pattern of zeros given 
by the continuum scaling limit of the lattice model. 

7 Conclusions 

The most important result we obtained in this paper is the state-by-state correspondence between 
the A4 ABF lattice model and the integrable QFT defined as the boundary (^1^3 perturbation of 
the c = 7/10 minimal unitary CFT describing the universality class of TIM. 

Such a result can be added to the previous corpus of known information and results from both the 
A4 lattice model and the corresponding c = 7/10 CFT possibly perturbed by the boundary operator 

^Our notation for the eigenvalues of T differs from that of [16] in order to avoid confusion with those of the 
normaUzed transfer matrix t. 
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(fi^s- As an example we can mention that in [7] it was found that the Hmit of the lattice model 
exactly reproduces the Virasoro characters of the corresponding c = 7/10 minimal unitary CFT. It 
is tempting to suggest that the presence of common structures (functional equations, characters) 
implies a strong (algebraic) identification between the A4 lattice model at the scaling limit and the 
conformal theory. 

Let us remember that the crucial ingredient which allowed us to establish the correspondence 
of states is the infinite set of integrals of motion coming from integrability. By means of them, 
an unambiguous procedure to fix the state-by-state identification has been given, and supported 
by a clear numerical evidence, at both the UV and IR critical points (the sectors (1,2) and (1, 1) 
respectively) establishing a true lattice-conformal dictionary (see Section 5). 

Furthermore, we followed the behaviour of the IM along the boundary flow. Since it is known 
that the integrable structure of {l2n-i{0}n survives off criticality, we can see, using the description 
of the flow given in [2] and the results of the present paper, that there is a smooth flow of the states 
which allows to see which state in the sector (1, 2) flows to which state in the sector (1, 1). 

It is important to stress that our constructive method to build the lattice-conformal dictionary 
essentially works level-by-level with an unambiguous procedure, however at present we are not able 
to see whether it is possible to write general expressions or not. Such a situation seems to resemble 
that of the computation of the expressions for the operators l2n-i' they can be computed one after 
the other but general formulas are not available. 

Another important result regards the connection (see sect. 6) between the "T-operators" 
defined in [16] with the transfer matrix t given in sect. 3.1. 

Let us conclude with some remarks about possible developments and extensions of the present 
work. Firstly, the existence of a lattice-conformal dictionary seems to suggest the possibility to 
define a "lattice Virasoro algebra" at least within the lattice ABF models (interesting approaches to 
a lattice Virasoro algebra can be found in [21, 22]). At present such a construction is prevented by 
the fact that we are not able to write down general expressions for the state-by-statc correspondence. 
Secondly, we expect one can get a clear idea of what quasi-particle states are at and off criticality. 

Finally, one can think to generalize the procedure to all the cases where a TEA is known for 
excited states. 
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Appendix 



Table Al. List of common cigcnstatcs and eigenvalues of I^, J5 in the vacuum sector of minimal models. The eigenstates are orthogonal. 
The first column is the level degeneracy (l.d.) as given by each single monomial appearing in the character expression. 



l.d. eigenstate 



h 



1 |0) 

W i^-2|0) 

L_3|0) 

2q^ [(I + 2c + V37 + 22c + 4c2)L_4 + ^L'^--^ |0) 
((1 + 2c - V37 + 22c + 4c2)L_4 + BLia) |0) 



c(22+5c) 
2880 



(480+c)(22+5c) 
2880 



(2160+c)(22+5c) 
2880 



96000+16342c+5c^+5760V37+22c+4c^ 
2880 



96000+16342c+5c^-5760V37+22c+4c^ 
2880 



2g5 ((5c + V 5(48 + 24c + 5c^) )L_5 + 2QL.,L.,) |0) ^13600+34822c+5c^+4320V5(48+24c+5c^) 
((5C - V 5(48 + 24c + 5c^) )L_5 + 2QL^,L_,) |0) ^13600+34822c+5c^-^4320V5(48+24c+5c^) 



-c(14+3c)(68+7c) 
290304 



(68+7c) (26208+9058c-3c^ ) 
290304 



(68+7c)(341712+104314c-3c^) 
290304 



(68+7c)(1141056+260050c-3c^+120960V37+22c+4c^) 
290304 



(68+7c)(1141056+260050c-3c^-120960V37+22c+4c^) 
290304 



(68+7c)(^4359600+869386c-3c2+151200-y/5(48+24c+5c2)) 

290304 

(68+7c)(^4359600+869386c-3c2-151200-y/5(48+24c+5c2)j 
290304 



4g6 (a(A)L_6 + AL-4L-2 + 7(A)^-3 + ^-2) |0) 
i = l,...,4 



12A + 83 + fiic+^c 



576"^ 



50-5c+12/3f 



310805 7746725c 205809050^ 238025c3 ^ 50* 



2592 



72576 



72576 ^ 13824 



+(87403 _ 212741, _ 5^ _ _ci_)^, + 80(56 + + 120^^ 



We used: a{f3) = ^2/33-6(i3c+8)/?^+2(i0ey_ic-^ ^^^^ ^ '^'"i^^rsffrof'""'^ ^ A are the 4 solutions of the following quartic equation: 

216/3^ - 54(15c + 16)13^ + 18(38c2 - 39c - 428)p^ - (IGOc^ - 1509c2 - 10998c - 4096)/3 - 12(30c3 + TOc^ - 923c - 1526) = 



Table A2. List of common eigenstates and eigenvalues of Is, I5 in the (1,2) sector of TIM, A = 1/10. The eigenstates are orthogonal. The 
first column is the level degeneracy (l.d.) as given by each single monomial appearing in the character expression. 



l.d. eigenstate 



1 10 



121 

' 19200 



1359 
512000 



^-lIlTj) 



26279 
19200 



668799 
512000 



196679 
19200 



19048239 
512000 



2g3 (|l6L_3 + (3 + V649)L_2L_i) | ^) 
(l6L_3 + (3 - V649)L_2-L-i) 1 1^) 

3g4 (L_4 + aiL_3L-i + /?(aOL2_2)|^) 
1 = 1,. ..,3 



427559 _|_ 9V649 



19200 



20 



427559 9V649 



19200 



4309553 
134400 



20 



+ 



216 



7(16-7ai) 



78251679 , 297V649 
512000 64 

78251679 297v^649 



512000 



64 



9(-864903608464- 106010360777a, +173389286400ai ^) 
512000(-43504+17953ai ) 



i = h...A 



-16904148+3744000Ki +82844337, 3(-2047389436+66944117l7i +788480007,^+752339200wi+1010688007jK;,) 
19200(-12+77i) 512000(-12+77f) 



We used: /3(a) = '^fsjie^yj)^ are the 3 solutions of the following cubic equation: 334a^ - 14210^ + 1072q; + 960 = 0; 

, . 2(50-357-872+115k-407k) , , ^1 ^ 1 j.- r xi r 11 • x / / 12 \ 

77(7, k) = — i5(_i2+7-y) 1 (7j) i^i) are the 4 solutions or the following system (7^7^-^]: 

r 84400 - 12910 7 - 5932972 + 191287=^ + 103300 k - 1169957 k + 47920 7^ k - 74750 + 26000 7 = 
\ 200-3807 + 10872 + 530/«-647k- 325k2 = 
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